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ABSTRACT 

We test the effectiveness of photometric redshifts based upon galaxy spectral template fitting for X-ray luminous 
objects, using a sample of 65 sources detected by Chandra in the field of the Caltech Faint Galaxy Redshift 
Survey (CFGRS). We find that sources with quasar-dominated spectra (for which galaxy spectral templates are 
not appropriate) are easily identified, and that photometric redshifts are robust for the rest of the sources in our 
sample. Specifically, for the 59 sources that are not quasar-dominated at optical wavelengths, we find that the 
photometric redshift estimates have scatter comparable to the field galaxy population in this region. There is no 
evidence for a trend of increasing dispersion with X-ray luminosity over the range Lx = 10 39 -5 x 10 43 erg s -1 , 
nor is there a trend with the ratio of X-ray to optical flux, fx/ fR- The practical implication of this work is that 
photometric redshifts should be robust for the majority (~90%) of the X-ray sources down to fx ~ 10~ 16 erg s" 1 
cm -2 that have optical counterparts brighter than R « 24. Furthermore, the same photometry can be easily used to 
identify the sources for which the photometric redshifts are likely to fail. Photometric redshift estimation can thus 
be utilized as an efficient tool in analyzing the statistical properties of upcoming large Chandra and XMM-Newton 
data sets and identifying interesting subsamples for further study. 



1. INTRODUCTION 

Photometric redshift estimation is a powerful tool in extra- 
galactic astronomy, and significant effort has been expended in 
recent years developing robust redshift estimation techniques. 
The most widely employed approach is to use spectral tem- 
plates (either empirical or from stellar synthesis models) to de- 
rive optimal fits to the observed galaxy colors (e.g. Lanzetta, 
Yahil, & Fernandez-Soto 1996; Fernandez-Soto, Lanzetta, & 
Yahil 1999; Benitez 2000; Furusawa et al. 2000). This ap- 
proach has been highly successful for normal galaxies, achiev- 
ing results as good as <j z = 0.06(1 +z) for the Hubble Deep Field 
(HDF; Fernandez-Soto^ Lanzetta, & Yahil 1999; Benitez 2000; 
Furusawa et al. 2000). Extending upon this work, Budavari et 
al. (2001) and Richards et al. (2001) have recently found that 
photometric redshifts accurate to within Az = 0.2 can be ob- 
tained for 70% of quasars in the Sloan Digital Sky Survey - 
even at z < 2.2, where quasar spectra lack a strong continuum 
break in the Sloan ugriz filters. 

In this work we focus on X-ray luminous objects, asking 
whether a self-consistent technique can be devised to obtain 
photometric redshifts for X-ray selected samples. Recent work 
on deep Chandra fields indicates that the resolved X-ray back- 
ground is comprised of a variety of objects, including early- 
type ellipticals, starburst galaxies, obscured active galactic nu- 
clei (AGN), narrow- and broad-line AGN, and quasars (see e.g. 
Mushotzky et al. 2000; Barger et al. 2001 ; Hornschemeier et al. 
2001; Tozzi et al. 2001; Stern et al. 2002; Barger et al. 2002). 
Given the disparate nature of these sources, the effectiveness of 
photometric redshifts is a priori unclear. Recent spectroscopic 
observations in the Chandra Deep Field-North (CDF-N) high- 
light this concern. Barger et al. (2002) find that, out of a sample 
of 182 hard sources (2-8 keV detections) with spectroscopic 



follow-up, approximately half show signatures of AGN activity 
in their spectra. 

As argued by Barger et al. (2002), one can reasonably expect 
that traditional galaxy spectral template fitting should be reli- 
able for the 50% of sources that lack an AGN spectral signa- 
ture. Similarly, it is likely that quasar spectral templates can be 
utilized for AGN-dominated sources (although this expectation 
has not yet been verified). It is unclear though whether either 
approach is robust for the intermediate sources in which both 
the AGN and the galaxy stellar population contribute signifi- 
cantly to the optical spectrum. Do photometric redshifts based 
upon galaxy spectral template fitting gradually degrade with 
increasing fractional luminosity contribution from the AGN, 
or do they remain robust until the AGN contribution reaches 
some critical level? The Hubble Deep Field, which is the 
canonical field for testing photometric redshifts, offers little in- 
sight. There are only six X-ray luminous sources with spectro- 
scopic redshifts in the HDF - the majority of which are low- 
luminosity. Photometric redshift comparisons in the HDF, such 
as the blind check of Cohen et al. (2000), thus cannot test the 
reliability of photometric redshift estimators in the interesting 
regime of properties. 1 

There are strong practical motivations for using an expanded 
sample of X-ray selected sources to assess the reliability of 
photometric redshifts for these objects. Foremost, an efficient 
means of redshift estimation is required to maximize the re- 
turn from large area Chandra and XMM-Newton surveys, such 
as the Lockman Hole survey (Chandra, PI: Barger), the Chan- 
dra Multiwavelength Project (ChaMP Wilkes et al. 2000), the 
XMM-LSS survey (Pierre et al. 2001), and the upcoming Chan- 
dra survey in the NOAO Deep- Wide Field (9 sq. degrees, Pis: 
Jones & Murray). A number of the issues that these surveys 
aim to address, such as evolution in the AGN X-ray luminosity 



1 One of these sources is an outlier in the blind photometric redshift tests of Cohen et al. (2000); however, it has also been argued by Fernandez-Soto et al. (2001) that 
the spectroscopic redshift for this source is incorrect. Two of the other objects are listed as faint X-ray sources in Hornschemeier et al. (2001), but do not appear in 
the «1 Ms Chandra catalogs of Brandt et al. (2001). 
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FlG. 1. — Optical R-band magnitude from Barger et al. (2002) plotted 
against the full band 0.5-8 keV X-ray flux from Brandt et al. (2001). The 
crosses correspond to the seven sources in our sample with colors that are best 
fit by a quasar spectral template. The filled circles correspond to the other 
sources in our sample, while the open circles denote the rest of the extra- 
galactic sources in the Barger et al. catalog that lie within the CFGRS region. 
Arrows denote sources that only have lower (upper) limits on their R-band 
magnitudes (X-ray fluxes). The dashed lines denote lines of constant fx /fit, 
as given by Equation 3 in Hornschemeier et al. (2001). 



function and evolution in the correlation function of faint X-ray 
sources, are statistical in nature and do not require the accu- 
racy of spectroscopic redshifts. For example, the survey in the 
NOAO Deep Wide Field aims to study the X-ray spatial correla- 
tion function using an expected ^2200 sources (C. Jones 2002, 
private communication). This project has no associated large 
spectroscopic program and will depend upon the use of pho- 
tometric redshifts to achieve its key scientific goal. If one can 
attain a level of accuracy for these photometric redshifts com- 
parable to what is achieved for "normal" galaxies, photometric 
redshifts will be a valuable tool. If not, then the limitations of 
photometric redshift techniques must be established. 

In this paper we assess the validity of photometric redshifts 
in the Caltech Faint Galaxy Redshift Survey field (CFGRS, 
Cohen et al. 2000; Hogg et al. 2000) using the Brandt et al. 
(2001) catalog of X-ray sources in the CDF-N. We first com- 
pare the scatter in photometric redshift estimates for X-ray se- 
lected sources with a control sample of quiescent galaxies, and 
quantify the fraction of X-ray sources with bright optical coun- 
terparts for which traditional galaxy spectral template fitting is 
valid. Next, we search for trends between the redshift residuals 
and other physical quantities, testing whether there is a range 
in X-ray luminosity (Lx) or flux ratio (fx/ fa) over which pho- 
tometric redshifts degrade. Finally, we ask whether the objects 
for which photometric redshifts are likely to fail (e.g. quasars) 
can be easily identified using the optical photometry. The data 
are described in §2 and the results are presented in §3. In §4 we 
summarize our work and discuss prospects for future studies 
with larger and fainter samples. We assume Oo = 0.3, = 0.7, 
and H = lOOh km s" 1 Mpc" 1 . 
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FIG. 2. — Template quasar spectra used in the photometric redshift code. 
The dashed line separates the SDSS composite spectra (short wavelength), 
which is corrected for stellar contamination, from the PDS 456 spectra. Simp- 
son et al. (1999) obtained /- and AT-band spectra for PDS 456; between these 
two bands we interpolate. 



2. DATA 

We focus on the CFGRS field due to the multiwavelength ob- 
servations and extensive spectroscopy available for this region. 
We include in our sample X-ray sources detected by Brandt et 
al. (2001) that lie within the CFGRS, have published spectro- 
scopic redshifts, and have unambiguous optical counterparts. 
The Brandt et al. (2001) source list is derived from the 975.3 
ks Chandra image of the Chandra Deep Field North (CDF-N), 
with minimum detectable fluxes near the aim point correspond- 
ing to f x « 1.9 x 10" 16 ergs s" 1 cm" 2 in the hard band (2-8 keV) 
and f x w 2.9 x 10" 17 ergs s" 1 cm" 2 in the soft band (0.5-2 keV). 
Spectroscopic redshifts (Cohen et al. 2000; Hornschemeier et 
al. 2001; Barger et al. 2002) and U„, G, TZ, K s photometry 
(Hogg et al. 2000) exist for 66 of these sources, which are listed 
in Table 1. For the majority of these objects the spectra are also 
published in Figure 6 of Barger et al. (2002). 

Of the 66 sources, only four are in the HDF proper and have 
previous published photometric redshifts using 7-color HST 
photometry. For one galaxy, HDF 36569_1302, Fernandez - 
Soto et al. (2001) dispute the spectroscopic redshift, arguing 
that the published value is attributable to the halo of a nearby, 
brighter galaxy. Due to this uncertainty, we exclude HDF 
36569_1302 from our analysis, yielding a final sample of 65 
sources. Figure 1 shows the distribution of optical and X-ray 
properties in our subsample, compared to properties of the rest 
of the extragalactic sources from Brandt et al. (2001) that lie 
within the CFGRS region. Our subsample includes 80% of 
the extragalactic sources with optical counterparts brighter than 
R = 24.5 (90% with counterparts brighter than R = 24). Since 
spectroscopic redshifts are lacking for sources with R > 24.5 
we are unable to study this subset of the population, which in- 
cludes the majority of faint X-ray sources with high f x / /r. We 
emphasize that the results presented in this paper are not neces- 
sarily applicable to these sources. 
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FIG. 3. — (a) Photometric redshift vs. spectroscopic redshift. The crosses are the sources that are AGN-dominated at optical wavelengths (§3), while the solid 
circles are the other sources in our sample. The dashed lines correspond to the lcr dispersion a, = 0.14(1 +z) measured for the Cohen et al. (2000) sample at z < 1.1. 
For the AGN-dominated sources, the redshift is based upon the best fit to the quasar spectral template discussed in §3. (b) Same as a, but with a more restricted 
redshift range (z < 1.7) and with the AGN-dominated sources excluded. 



3. PHOTOMETRIC REDSHIFTS 
3.1. Redshift Estimation Technique 

We utilize the Bayesian photometric redshift code of Ben- 
itez (2000), BPZ, to derive photometric redshifts for our sam- 
ple. This code employs spectral template fitting, but also takes 
into account the /-band magnitude of the galaxy. Specifically, 
the code utilizes the F814 magnitude distribution of the Hub- 
ble field as a Bayesian prior probability distribution, which en- 
ables superior discrimination when color degeneracies exist be- 
tween low- and high-redshift galaxies. For the current analysis 
we have only four filters, so the Bayesian magnitude prior is 
necessary to break these color degeneracies and minimize the 
number of spurious outliers. Our results should be valid for any 
code that utilizes spectral template fitting though, and to verify 
this expectation we also derive photometric redshifts using an- 
other publicly available code, Hyperz (Bolzonella, Miralles, & 
Pello 2000). While BPZ is more robust to outliers for our data 
set, we confirm with Hyperz that our results are insensitive to 
which photometric code is used. 

The standard six spectral templates described in Benitez 
(2000) are used for this analysis. Four are based upon spec- 
tral energy distributions (SED's) in Coleman, Wu, & Weedman 
(1980) (E/S0, Sbc, Scd, and Irr), and two are derived from spec- 
tra for starburst galaxies in Kinney et al. (1996). One challenge 
for the current analysis is that we lack /-band data for compar- 
ison with the HDF magnitude priors. Fortunately, the code is 
not strongly sensitive to errors in the /-band magnitude used for 
the prior, so we make the coarse assumption that TZ-F%\A ~ 1 
(comparable to an elliptical galaxy at z — 0.5) for all galaxies 
in the sample. We also test the code with IZ-FRIA = 0, finding 
that the choice of fiducial color yields a negligible difference in 
estimated redshifts. 



We employ a two-part approach in deriving photometric red- 
shifts. A necessary first step is to identify objects with quasar- 
dominated spectra at optical wavelengths, since photometric 
redshifts based upon galaxy spectral templates are likely to fail 
for these sources. If these objects can be identified, then pho- 
tometric redshifts can potentially be derived for them by other 
means (e.g. see Budavari et al. 2001; Richards et al. 2001). To 
this end, we construct a quasar template using the optical SDSS 
composite spectrum of Vanden Berk et al. (2001), coupled with 
infrared spectra of PDS 456 from Simpson et al. (1999). For 
the SDSS composite, we apply a correction for stellar contami- 
nation, assuming that the contamination increases linearly from 
10% at 6000A to 30% at 8500A, which is based upon estimates 
in Vanden Berk et al. (2001). The resulting template is shown 
in Figure 2. 

We initially include this quasar spectral template in the 
SED library and run the photometric redshift code without the 
Bayesian prior. Six of the 65 sources in our sample, including 
the two highest redshift sources, 2 are best fit by the quasar spec- 
trum. One of these sources has previously been identified as a 
quasar (Liu et al. 1999; Hornschemeier et al. 2001); a second 
was identified by Hornschemeier et al. (2001) as a BL AGN. 
We exclude these six sources from the subsequent analysis, but 
record the photometric redshift estimates and include them in 
subsequent Figures. For the remaining 59 X-ray sources, all of 
which are at z < 1 .3, we perform a second iteration of the red- 
shift estimation in which we exclude the quasar spectrum and 
include the Bayesian magnitude priors. The results are shown 
in Figures 3a, b. We derive a scatter a z = 0.10(1 +z) for these 
sources. 



2 This is a selection effect - if a source at z — 2 is not a quasar, then it will be optically faint and likely not have been targeted for spectroscopy. 
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FIG. 4. — Comparison of photometric redshifts obtained with Hyperz and 
BPZ for the X-ray luminous sources. We only plot objects that are neither 
quasar-dominated nor outliers with Hyperz (as defined in the text). 

3.2. Comparison with Normal Galaxy Sample 

Cohen et al. (2000) provide spectroscopic redshifts for a total 
of 529 galaxies in the Hogg et al. (2000) photometric catalog, 
enabling construction of an X-ray quiet comparison sample. 
Excluding the Chandra sources and galaxies at z > 1.3, the red- 
shift dispersion for this catalog is a z = 0.18(1 + z). This scatter 
is dominated by a handful of outliers - if we employ er-clipping 
with rejection of outliers beyond 3.5a, then 1 1 out of 461 galax- 
ies are rejected and the dispersion reduces to a z = 0.14(1 +z). 
This value is larger than the scatter for the HDF (cr, < 0. 10(1 +z) 
for 7 filters; e.g. Fernandez-Soto, Lanzetta, & Yahil 1999; Ben- 
itez 2000; Bolzonella, Miralles, & Pello 2000; Furusawa et al. 
2000), but not unreasonable given the larger photometric un- 
certainties and smaller number of passbands in this analysis. 
Comparing with the X-ray subsample, we conclude that, with 
the exception of clearly quasar-dominated sources, the disper- 
sion in photometric redshifts is no greater for X-ray sources 
than for the field galaxy population. 

To verify these results, we also derive photometric redshifts 
using Hyperz. The Hyperz code was utilized by Barger et al. 
(2002), who concluded that the lack of priors coupled with their 
lack of {/-band data significantly hindered the reliability of 
their photo-z's. Our comparison of BPZ and Hyperz confirms 
that the inclusion of a magnitude prior significantly improves 
the robustness, particularly with a limited number of passbands. 
With Hyperz, 12% of the sources in the quiescent sample - and 
a comparable fraction of the X-ray luminous sources - are out- 
liers with Az > 1 . In contrast, only 2% of sources in the qui- 
escent sample are outliers with BPZ. Nonetheless, if these out- 
liers are excluded, then Hyperz and BPZ typically yield consis- 
tent photometric redshifts for individual X-ray selected sources 
(Figure 4). 3 Furthermore, Hyperz yields dispersions compara- 
ble to BPZ, with values a z = 0.14(1 +z) for the full field sample 
and a z = 0.13(1 +z) for the X-ray sources. Both codes thus 

3 Hyperz does appear to underestimate the redshift for sources at z ~ 0.7 though. 
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FIG. 5. — Photometric redshift error as a function of X-ray luminosity. 
The symbols are the same as in Figure 3. The dashed lines correspond to 
a, = 0.14(1 +z). The vertical axis is set to match the dynamic range of the 
non-quasar-dominated sources. Consequently, two highly discrepant quasar- 
dominated sources with dz/(l +z) ~ 1.5 and L x ~ 4 X 10 42 erg s _1 cnr 2 are 
not visible in this plot. 
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FIG. 6. — Photometric redshift error as a function of fx/ /r- The symbols 
are the same as in Figure 3, and the dashed lines correspond to cr z = 0. 14(1 +z). 
As in Figure 5, two quasar-dominated sources with dz/(l+z) ~ 1.5 and 
fx I fit ~ 2 - 4 are not visible in this plot. 

indicate that the dispersion is comparable for the X-ray lumi- 
nous and quiescent samples. 

3.3. Dependence upon X-ray Properties 

The results of the preceding sections indicate that photo- 
metric redshifts are robust for ^90% of X-ray sources with 
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R < 24.5 down to fx ~ 10~ 16 erg s -1 cm" 2 , including some 
sources with a moderate AGN contribution. One naively ex- 
pects increased scatter in the photometric redshifts with increas- 
ing optical contribution from the AGN. A direct test of this ex- 
pectation would require fitting a linear superposition of normal 
galaxy and quasar spectral templates to obtain the fractional 
contribution, which is beyond the scope of the current analysis. 
Instead, we ask whether the probability of failure is correlated 
with the X-ray properties of the source - specifically the X-ray 
luminosity and the ratio of X-ray to optical flux {fx/ fid- The 
latter quantity is the best analog measure of the relative AGN 
and starlight contributions, but with the caveat that there exists a 
factor of ten range in fx /fit within which sources can be either 
AGN or starlight dominated. 

In Figure 5 we plot dz/(l +z) against L x , where d z = z p hot - 
z S pec, to assess whether the redshift estimates degrade above 
some threshold luminosity. Aside from objects that are best 
fit by the quasar spectral template, there is no indication of 
degradation up to Lx = 5 x 10 43 /z 2 erg s" 1 , the luminosity of 
the brightest non-quasar-dominated sources in our sample. We 
also find no statistically significant trend with fx/ fit (Figure 
6); however, the probability of identifying a source as quasar- 
dominated increases significantly for sources with fx/ fit <^ 0.5. 
Nearly 30% of sources above this threshold (5 of 18) are best 
fit by the quasar template, as opposed to only one of 47 at lower 
fx /fit. These results imply, somewhat surprisingly that the 
photometric redshift estimates are fairly insensitive to the pres- 
ence of an AGN unless the AGN dominates the optical spec- 
trum. 

4. SUMMARY 

We use a sample of X-ray sources detected by Chandra in 
the Caltech Faint Galaxy Redshift Survey region to test the ro- 
bustness of photometric redshifts for these objects. For the 59 



out of 65 sources with colors that are not best fit by a quasar 
spectral template, we find no degradation in the accuracy of 
photometric redshifts as compared to the field galaxy popula- 
tion. We also find that the redshift residuals exhibit no trend as 
a function of X-ray luminosity or fx/ fit for these sources. We 
demonstrate that it is feasible to quickly identify objects whose 
spectra are quasar-dominated at optical wavelengths and derive 
robust photometric redshifts for the other ~ 90% of sources that 
have optical counterparts brighter that R s=s 24.5, which includes 
roughly two thirds of all sources with fx > 10~ 16 erg s" 1 cm" 2 
(Barger et al. 2002). Consequently, photometric redshift esti- 
mation should be a valuable tool for upcoming large samples of 
optically bright, X-ray selected objects. This paper is a first step 
towards quantifying the robustness of photometric redshifts for 
systems in the transition regime between quiescent galaxies and 
quasars. The next required step is extension to fainter optical 
magnitudes. Given that the majority of fainter sources have 
high fx /fa, it is unclear whether photometric redshifts will be 
as effective for this population, but the motivation for testing 
their robustness is strong. There is speculation that these opti- 
cally faint, high fx / /r sources possibly include obscured AGN 
at z > 5 (e.g. Stern et al. 2002), and photometric redshifts may 
be the only viable approach for studying this population. 
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Comments* 5 



QSO (Liu et al. 1999) 



z«=0.10 

z FS = 0.85 
zfs = 0.38 



Excluded; zfs = 1 -27 



"Identification number from Brandt et al. (2001). 
b Identification number from Hornschemeier et al. (2001). 

c Full band flux (0.5-8 keV) in units of 10~ 15 erg s~' cnT 2 from Brandt et al. (2001). 

d Spectroscopic redshifts are taken from Barger et al. (2002). The one exception is source 220, for which 
the redshift is taken from Cohen et al. (2000). 

e 'Q' denotes that the photometry for this source was best fit by a quasar spectral template. 

f 'Y' denotes that this source was found to have broad emission lines in the spectra by Barger et al. (2002). 



